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Abstract 

In adsorption study, to describe sorption process and evaluation of best-fitting isotherm model is a key analysis to 
investigate the theoretical hypothesis. Hence, numerous statistically analysis have been extensively used to estimate 
validity of the experimental equilibrium adsorption values with the predicted equilibrium values. Several statistical 
error analysis were carried out. In the present study, the following statistical analysis were carried out to evaluate 
the adsorption isotherm model fitness, like the Pearson correlation, the coefficient of determination and the 
Chi-square test, have been used. The ANOVA test was carried out for evaluating significance of various error 
functions and also coefficient of dispersion were evaluated for linearised and non-linearised models. The adsorption 
of phenol onto natural soil (Local name Kalathur soil) was carried out, in batch mode at 30 ±20 C. For estimating 
the isotherm parameters, to get a holistic view of the analysis the models were compared between linear and 
non-linear isotherm models. The result reveled that, among above mentioned error functions and statistical 
functions were designed to determine the best fitting isotherm. 
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Introduction 

Phenols are primarily present in oil refinery, coal, paper, 
textile, synthetic rubber and pharmaceutical wastewaters 
and are an problem of a serious environmental due to their 
high toxicity and potential accretion in the environment. 
There are numerous methods for removing phenols from 
the wastewaters, such as reverse osmosis, adsorption, 
bio-degradation, chemical oxidation and solvent extraction. 
Adsorption is an efficient treatment system for removing 
phenols from wastewaters. 

Adsorption isotherms experimental data is necessary 
for the design of sorbate-sorbent. To remove phenols 
from wastewaters a definite sorbate-sorbent system has 
required to optimize the design and also experimental 
equilibrium data is important to establish the most 
appropriate correlation. Kumar and Sivanesan has been 
studied comparison of linear and non-linear methods 
for removal of safranin onto rice husk. Subramanyam 
and Das has been studied comparison of linear and 
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non-linear isotherm models, adsorption of phenol 
onto natural soil (Kalathur soil). The various adsorption 
isotherm equation to remove pollutants from wastewaters 
as have been used to study the nature of adsorption, 
with the fundamental idea of optimization of the design 
parameters. Many most commonly used isotherms models 
appearing in the adsorption literature including, Langmuir, 
Freundlich and Redlich-Peterson isotherm models [1,2]. 
All these isotherm models derived based on the theoretical 
assumption, to measure the goodness of fit in the literature 
there are number of statistical error deviation functions 
such as the Marquardts percent standard error deviation 
(MPSED), the correlation coefficient (r^), the sum of the 
squares of the errors (SSE), the hybrid fractional error 
function (HYBRID), the average relative error deviation 
(ARED) and the residual analysis (RESID) [3]. However, 
the very approach of linearization of the nonlinear models 
necessarily yields rationalization of specific variables, 
which may have significant bearings on the adsorption 
process itself. 

The main objective of this study was to explore, the 
applicability of the statistical methods in determining 
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the best fitting isotherm models. The statistical tools used 
in the present study were Pearson correlation coefficient 
(r), the coefficient of determination (r^), the Chi square test 
(X^) and ANOVA test. The single component sorption 
study was carried out based on experiment and computed 
phenol-sorption by the soil selected (namely, Kalathur soil) 
on phenol from its aqueous solution to estimate the 
effect of linearization on the accuracy of the model 
(on comparison with experimental values). 

Material and methods 

The soil (namely Kalathur soil) sample was collected 
from Thanjavur districts, Tamil Nadu (India). The 
samples, thus obtained, was washed thoroughly with 
distilled water and dried for 2 hours, at 105°C in an 
electric oven, followed by crushing and sieving (100-635 
SIEVE NO ASTM Ell-87), to obtain the uniform size of 
particles. The final sample, after passing through the sieve, 
was dried, desiccated and preserved in air-tight chamber 
for subsequent analysis and experiments. Soil texture ana- 
lysis was carried out to find the percentage of sand (25%), 
silt (20%) and clay (55%) present in the soil. According to 
United States Department of Agriculture (USDA) texture 
triangle Kalathur soil (Kr) was classified as clayey soil. 

Phenol (CeHsOH) of analytical reagent (AR) grade 
supplied by Ranbaxi Laboratories Ltd., India, was used 
for the preparation of synthetic adsorbate of concentra- 
tion 100 mg/1. The required quantity of phenol was 
accurately weighed and dissolved in distilled water and 
make up to one liter. Fresh stock solution was prepared 
every day and stored in a brown color glass bottles to 
prevent photo-oxidation. 

To study the equilibrium, batch experiments were con- 
ducted at room temperature (30 ± 2°C) for an adsorption 
period of 24 hours. The effect of adsorbent dosage on 
the uptake of phenol on to the soil (namely Kalathur 
soil) was studied at different adsorbent doses (50 to 
1000 mg/100 ml) for the concentration of 100 mg/1. The 
percentage phenol removal and equilibrium adsorption 
uptake, qe (mg/g), was computed by making use of the 
equation as given below: 

Percent removal = 



lOO(Co-Q) 
Co 

Adsorbed amount (mg/g) 

_{Co-Ce)V 



w 



(1) 



(2) 



Where Cq is the initial phenol concentration (mg/1), 
Ce the equilibrium phenol concentration (mg/1), V the 
volume of phenol solution (1) and w is the soil mass of 
the adsorbent (g). 



Estimation of best-fitting isotherm model 
Error functions 

Average relative error deviation (ARED) is to minimize 
the distribution in fractional error over the entire range 
of concentration studied [3]. 



ARED 



^e.cal Qe. 



Ce.exp 



X 100 



(3) 



This error function has a major drawback, inspite of that 
most of the researchers prefer using this error function. At 
higher end of the liquid state concentration rage, the 
calculated isotherm parameters obtained from such error 
function will yield a better fit. This has been resulted 
because of the magnitude of the errors and therefore the 
error function will increase as concentration increases. 



SSE = Y,{Qe,carQ< 



e.exp ) 



(4) 



The hybrid fractional error function (HYBRID): To 
improve the sum of the squares of the errors at lower 
levels of liquid-phase concentration, this error function 
was developed. In this task, each the sum of the squares 
of the error values was divided by the theoretical adsorbent 
phase concentration value. 



HYBRID 



N-P 



E 



ce.exp 



X 100 



(5) 



The Marquardts percent standard deviation (MPSED): 
This error function distribution follows the geometric mean 
error which allows for the number of degrees of freedom of 
the system. 



MPSED - 



\ 



^ Qexp Qcal ) 



N-P 



-XlOO 



(6) 



The sum of the absolute errors (EABS): It is similar to 
SSE and provides a better fit at higher concentration for 
the isotherm parameters. 

p 

EABS=J2\Qe,exp-Qe,cal\ (7) 

i=l 

Statistical functions 

Correlation coefficient of Pearson (r): It is a sampling 
index, shows the degree linearity of between two dependent 
data series. The degree of linearity varies from -1 to 1. 



^E^^-(E^)' ^E^^-(Ei') 



(8) 
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Table 1 Linear and non-linear isotherm equation 



Isotherm 



Non-Linear 
equation 



Linear equation 



Langmuir-l 




Q _ Q , ^ 

qe qm 


Langmuir-2 


_ qr. bCe 
He — i+bCe 


i = (b^) i + i 


Langmuir-3 






Langmuir-4 




J = bq^-bQe 


Freundlich 


qe = KpCy^ 


log(qe)= log(KF)+l|og(Ce) 


Redlich-Peterson 


- [l+a,p(Q)^] 


In(^-l) =pin(Ce)+ In(aRp) 



Coefficient of determination (r^): It explains the regression 
line with percentage of variability in the dependent data 
series variable. The percentage degree varies from 0 to 1. 



^{XY)<^{YY) 



(9) 



Where Sxy is the sum of squares of X and Y, Sxx is 
the sum of squares of X and Syy is the sum of squares 
of Y. In addition to above mentioned error and statistical 
functions Chi-square test, was also examined to predict 
best-fitting isotherm models. 

The ANOVA test (two factors without replication) 
was carried out for evaluating significance of various 
error-functions and four coefficients of dispersion 
(namely, Coefficient of Range, Coefficient of Quartile 
Deviation, Coefficient of Mean Deviation, Coefficient 

Table 2 Linearised isotherm parameters 



Isotherm 



Kalathur soil 



Langmuir-l 

qm 

b 

Langmuir-2 

qm 

b 

Langmuir-3 

qm 

b 

Langmuir-4 

qm 

b 

Freundlicli 

KF 

n 

Redlicli-Peterson 
Krp 

P 

a 



52.63 
0.044 

41.67 
0.068 

46.52 
0.055 

49.37 
0.049 

5.801 
2.747 

30.345 

0.987 

0.789 



Table 3 Non-linearised isotherm parameters 

Non-linearized isotherm 



Kalathur soil 



Langmuir 




qm 


51.83 


b 


0.04333 


R2 


0.9952 


Freundlicli 




KF 


5.635 


n 


2.175 


R2 


0.9953 


Redlicli-Peterson 




Krp 


2.351 


P 


0.9634 


a 


0.05369 


R2 


0.9953 



of Variation) were evaluated, separately for linearised and 
non-linearised models. In case of linearised models, only 
one type of Langmuir distribution (i.e., type-1) was 
considered because of low mean and lowest variance. A 
paired t-test was also carried out between the dispersion 
coefficients of linearised and non-linearised models to 
evaluate the t-statistics. 

Results and discussion 

In the present study, to find out the isotherm models 
(linear and non-linearized isotherm) that can describe 
with precision the experimental results of adsorption iso- 
therms compare the parameters that can be determined 
(linear and non-linearized isotherm) and also determine 
the theoretical adsorption isotherms. To remove phenol 
from liquid phase it is necessary to develop a relationship 
between a sorption-sorbate system and equilibrium data. 
Three most commonly used isotherms (viz. Langmuir, 
Freundlich isotherm and Redlich-Peterson equation) were 
studied. Table 1 shows isotherm models that are used in 
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Figure 1 Linearised iostherm models for phenol adsoprtion by 
Kalathur soil. 
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Figure 2 Non-Linearised iostherm models for phenol adsoprtion 
by Kalathur soil. 



the present work and their associated parameters are 
given Tables 2 and 3 for both linear isotherm and non- 
linear isotherm analysis. Figure 1 shows the fitting values 
of linear regression and Figure 2 shows fitting values of 
non-linear regression analysis. 

Incorporation of the ANOVA study for the error func- 
tions with regard to the isotherm models (for both linear 
as well as non-linear models) to estimate the significance 
of variance (refer to Tables 4 and 5). The model compatibil- 
ity (for both linear and non-linear models) was estimated 
for the error functions, with regard to not only the lowest 
value of error, but also the other coefficients of dispersion 
(refer to Table 4). 

Error analysis 
Linear isotherm 

The linearized four forms of Langmuir isotherm model 
were presented in Table 1 [4] and the Langmuir coefficients 
for four linearized Langmuir equation were obtained 
by plotting graphs between Ce/qe versus Ce (Type- I 
linearized equation), 1/qe versus 1/Ce (Type- II linearized 
equation), qe versus qe/Ce (Type- III linearized equation). 



and qe/Ce versus qe (Type- IV linearized equation). 
Table 2 shows the calculated parameters of the four 
linearized Langmuir isotherm model. A graph was 
drawn between experimental and observed data as 
shown in Figure 2. From Table 2, it can be inferred that, 
different linear Langmuir equations show different 
Langmuir constants, as indicated by variation in errors, 
specific to the corresponding mode of linearization [5]. In 
the case of Kalathur soil, on comparison of the four linear- 
ized Langmuir equations, it is observed that the Type- I 
linearized Langmuir equation showed higher value of 
correlation coefficient (r2 = 0.994) than that of the other 
three linearized equations (Type- II to IV) as shown in 
Table 2. The adsorption capacity of Kalathur soil was found 
to be 52.63 mg/g for Type -I linearized Langmuir and that 
of Type- II, III and IV are 41.67 mg/g, 46.52 mg/g and 
49.37 mg/g, respectively. Thus, during linearization, errors 
in the computation of parameters may be responsible for 
the variation in adsorption capacity q^' and adsorption 
constant 'b' (from Type-I linearized isotherm form to 
Type-IV linearized isotherm). In other words, the trans- 
formation of non-linear isotherm model to linear isotherm 
models seems to implicitly alter the error functions as well 
as the error variance and normality assumptions of the 
least-squares methods [4,6,7]. As suggested in the lower 
correlation coefficient values, it will be inappropriate to use 
this type of linearization. 

In order to verify the validity of the linearized and 
non-linearised isotherm models as well as the best-fitting 
isotherm model, six common statistical error methods 
were employed to calculate the error divergence between 
observed and predicted sorbate-sorbent system data. It is 
clear from the Table 6 that the linearized Langmuir models 
(Type-I, Type-II, Type-Ill and Type-IV), Freundlich and 
Redlich-Peterson isotherms, Langmuir isotherm Type-1 
shows higher r^ value and low error values (i.e., ARED, 
HYBRID and EABS). Which indicates among all the 



Table 4 Variation of mean & dispersion coefficients between linearised and Non-linear models 





Errors- > statistics 


ARED 


MPSED 


HYBRID 


SSE 


RESID 


EABS 


Linearised model 
















mean 




22.0 


27.9 


22.0 


437.0 


10.2 


35.4 


Dispersion 


Coefficient of Range 


0.9 


0.9 


0.9 


1.0 


2.6 


0.9 




Coefficient of Quartile Deviation 


0.5 


0.5 


0.5 


0.9 


-3.6 


0.7 




Coefficient of Mean Deviation 


1.0 


1.0 


1.0 


1.3 


3.5 


0.9 




Coefficient of Variation 


1.5 


1.5 


1.5 


1.9 


5.5 


1.2 


Non-linear model 
















mean 




4.9 


7.5 


4.9 


11.2 


0.1 


7.3 


Dispersion 


Coefficient of Range 


0.3 


0.1 


0.3 


0.7 


-8.4 


0.5 




Coefficient of Quartile Deviation 


0.2 


0.1 


0.2 


0.5 


1.8 


0.3 




Coefficient of Mean Deviation 


0.3 


0.1 


0.3 


0.8 


5.9 


0.5 




Coefficient of Variation 


0.4 


0.2 


0.4 


1.0 


7.7 


0.6 
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Table 5 ANNOVA 



Source of variation df Sum of squares Mean sum of square F-statistics P-value 

Linearised Non-linear Linearised Non-linear Linearised Non-linear Linearised Non-linear 

Among isotherms 2 600394.12 134.21 300197.06 67.11 1.74 3.73 0.23 0.06 

Between error Functions 5 1627971.12 204.89 325594.22 40.98 1.88 2.28 0.18 0.13 

Residual Error 10 1727359.62 179.72 172735.96 17.97 



linearised isotherm form of the Langmuir isotherm 
type-I is able to describe equilibrium data and indicates 
the best linearised isotherm model 

Redlich-Peterson isotherm (using three-parametric 
modeling equation: Table 2) was plotted using experimental 
data between ln(KrpCe/qe -1) and In(Ce). Redlich-Peterson 
isotherm contains three unknown variables (viz. Krp, a 
and |3), it is not possible to obtain three unknown variables 
using linearising Redlich-Peterson isotherm. Thus, the three 
unknown variables were obtained by minimization of 
the isotherm equation (and, thus, by maximization of 
the correlation coefficient) (Tables 6 and 7). In this 
case, the calculated parameters need not be unique 
(and could reflect the local optima) and hence the 
comparison of linearised and non-linear isotherms may 
not be relevant. In order to verify the model validity the 
correlation coefficient was lower, as well as the ARED, 
HYBRID, EABS value are very high (when compare other 
four linearized Langmuir and Freundlich). Therefore, the 
Redlich-Peterson linearised isotherm model fails to explain 
the sorbate-sorbent system of phenol onto Kalathur soil. 

Non-linear isotherm 

In the present study, for studying non-linear isotherm 
models. Graph Pad Prism versions 5.0 have been used 
for determining the non-linear coefficients. The determined 
coefficients were shown in Table 3. Correspondingly, a plot 
was drawn between Ce versus qe (Figure 2) using the 

Table 6 Error functions 

Error function type/isotherm model 

Linear 
Langmuir 
Type-! 
Type-!! 
Type-Ill 
Type-IV 
Freundlich 
Redlich-Peterson 
Non-Linear 
Langmuir 
Freundlich 
Redlich-Peterson 



experimental and predicted value by non-linear models. 
From Tables 6 and 7, it was observed that the correlation 
coefficient value is high (r ^ = 0.9953) and the low ARED, 
HYBRID and EABS values, thus it indicates that the 
models are able describe equilibrium data perfectly. 
Therefore, as far as the non-linear isotherm model is 
concerned the error remains constant. Hence, to use the 
correlation coefficient values for comparing the best-fitting 
non-linear isotherm models is befitting. 

Table 3 shows a non-linear Langmuir, Freundlich and 
Redlich-Peterson model parameter values and Figure 2 
shows a plot between Ce versus qe . The correlation 
coefficient value higher than that of linearized isotherm 
model. Table 6, shows error function value. But, the error 
values show the improved up on linearized isotherm. 

In the case of Langmuir and Redlich-Peterson model, 
the model shows high correlation coefficients value. The 
error functions ARED, HYBRID and EABS was found to 
be good for Langmuir isotherm model and error functions 
MPSED and SSE was found to be good for Redlich- 
Peterson isotherm models. It is clear from above results 
the Langmuir isotherm and Redlich-Peterson model were 
for better than Freundlich isotherm model. The models 
were able describe experimental data perfectly. Hence, it 
can be understood that, the Redlich-Peterson and 
Langmuir isotherms were the most suitable models for 
sorbate-sorbent system. A close correspondence was found 
to exist between Langmuir and Redlich-Peterson isotherm 



ARED 



MPSED 



HYBRID 



SSE 



RESID 



EABS 



5.156953 

8.543072 
6.432829 
5.219617 
19.68531 
86.76697 



7.148837 
9.621322 
7.416404 
6.884058 

23.10193 
113.0963 



5.156953 

8.543072 
6.432829 
5.219617 
19.68531 
86.76697 



7.578923 
56.61154 
16.64189 
7.380351 

481.0593 
2053.002 



4.022454 
-6.95591 
-1.62588 

0.566875 
-51.2595 
116.6 



7.469256 

1 7.98949 
1 1 .22502 
7.610797 
51.58342 
116.6088 



3.722937 

7.189632 
3.777602 



7.060828 
8.87458 
6.685821 



3.722937 

7.189632 
3.777602 



4.801368 
24.15513 
4.70022 



0.716535 
-0.91032 

0.541658 



4.474066 

12.61618 
4.739453 



The bold letters Indicates lowest values of ARED, MPSED, HYBRID, SSE, RESID and EABS. 
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Table 7 Statistical function 



Statistical function 
type/isotherm model 



Pearson 

correlation 

coefficient 



Determination 
coefficient 



Chi square 



Linear 

Langmuir 

Type-I 

Type-ll 

Type-Ill 

Type-IV 

Freundlich 

Redlich-Peterson 

Non-Linear 

Langmuir 

Freundlich 

Redlich-Peterson 



0.9419 
0.9149 
0.9287 
0.9357 
0.9759 
0.7951 

0.9485 
0.8955 
0.9453 



0.994 

0.951 
0.891 
0.891 
0.979 
0.98 

0.9952 
0.9953 
0.9953 



The bold letters indicates highest values of Pearson Correction 
determination Coefficient as well as lowest Chi Square. 



0.388526 

0.598918 

2.003552 

0.78796 

20.96659 

60.37647 

0.388526 
1.246057 
0.368326 
coefficient, 



models. Similar findings have also been reported by other 
researchers as well [3,4,8,9] . 

As can be seen fi:om Table 7, best fitting isotherm model 
was determined more appropriately by the Chi-square 
test. However, among the three isotherm models studied 
in this work, the phenol adsorption onto soil system was 
appropriately explained by Redlich-Peterson isotherm 
model only. Indeed the transformation of non-linear 
to linear models misrepresents the experimental error, 
which limits the validity of the error function and statis- 
tical tools. To avoid the errors discussed above, non-linear 
regression method can be more appropriately used. 

As per the ANOVA table, although the p-values for 
isotherms is higher in case of linearised model than in case 
of non-linear models, indicating the relatively higher signifi- 
cance level associated with non-linearised case than the 
linearised case, yet at 0.05 significant level the variation is 
less than the table values for F-distributions. In fact, the 
variability between error fiinctions is even less significant. 
These indicate that the selection of error fianctions and the 
isotherm models are fairly unbiased estimators for line- 
arised and nonlinear models (Tables 4 and 5). 

The overall mean of the errors of the linearised models 
is 93.5% higher than that of the non-linear models. To 
evaluate dispersion from the mean value, four coefficients 
of dispersion were studied, which indicated that 87.5% time 
the linearised models show higher dispersion coefficient. In 
fact, paired t-test results show that the confidence level for 
the non-linear models (as higher than linearised models) is 
more than 84%. 

Conclusions 

Thus based on the statistical studies it was found that 
variability in both linearised and non-linear cases are not 



significant at 0.05 significant level (both among isotherms 
and among error fianctions) confirming the error-estimators 
& isotherm models used as fairly unbiased (yet non-linear 
models do have relatively higher significance compared to 
linearised models, though). The overall mean of the error 
fianctions of linearised models is significantly higher (and so 
also most of the dispersion coefficients), compared to their 
non-linear counterparts, indicating non-linear modeling to 
be much better representation of experimental results than 
the linearised ones. 
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